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The three pentaquark candidates + (153O), 5(1862) and 0^(3100) have been stud- 
ied in ep collisions at a centre-of-mass energy ^fs = 300 — 318 GeV using the full 
luminosity of the HERA-I data. Searches for narrow baryonic states in the decay 
channels K^p, K+p, S — n^, S+tt^ and D* ± p = F are reported. The results support 
the existence of a narrow resonance decaying into K^p and Kgp, consistent with 
the 0(1530) state. No signals are seen in the K+p, H - ^, E+7T and 
channels. 



1. Introduction and Experiment 

The HERA e-p collider accelerates electrons (or positrons) and protons 
to energies of E e = 27.5 GeV and E p = 920 GeV (820 GeV until 1997), 
respectively. The two collider experiments, HI and ZEUS, are located at 
two collision points along the circulating beams. The incoming e^- interacts 
with the proton by first radiating a virtual photon. The photon is either 
quasi-real with Q 2 < 1 GeV 2 and Q 2 med%an ~ 3 • 10~ 4 GeV 2 , where Q 2 is 
the negative squared four-momentum transferred between the electron and 
proton, or highly virtual (Q 2 > 1 GeV 2 ). In the former case no scattered 
electron is visible and this is the photoproduction (PHP) regime. In the 
latter case the scattered electron is measured in the main detector and this is 
the deep inelastic scattering (DIS) regime. The analysis was performed with 
ZEUS data taken between 1995 - 2000 ("HERA-I"), corresponding to an 
integrated luminosity of « 120 pb _1 . Charged particles are tracked in the 
Central Tracking Detector (CTD) covering polar angles of 15° < 9 < 164°. 
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The energy loss of particles in the CTD, dE/dx, is estimated from the 
truncated mean of the anode-wire pulse heights, after removing the lowest 
10% and at least the highest 30% depending on the number of saturated 
hits. The dE/dx resolution for full-length tracks is about 9%. 

2. Evidence for the strange pentaquark ©+(1530) 

Fixed-target low-energy experiments saw a narrow exotic baryon with 
strangeness +1 around 1530 MeV decaying into K + n. It was attributed to 
the + = uudds pentaquark candidate predicted by Diakonov et al. 1 at 
the top of a 51/(3) spin 1/2 anti-decuplet of baryons. Narrow peaks were 
also seen at a similar mass in the final state K s p, which is not necessarily 
exotic. They were attributed to the 9 + as well. It is interesting to search 
for the + baryon in high-energy collider experiments. In particular it can 
be searched at the central rapidity region, which has little sensitivity to 
the proton remnant region. This region is dominated by parton fragmen- 
tation with no net baryon number, unlike low-energy experiments, where 
the pentaquark is mainly produced in the nucleon fragmentation region. 

The ZEUS search for the 6+(1530) 2 used the 1996 - 2000 HERA data 
(121 pb -1 ) and was performed in the DIS regime (Q 2 > 1 GcV 2 ). The 
search in the K°p mode was complicated due to a few unestablished reso- 
nances, such as £(1480) and £(1560), called "£ bumps" 3 . There are no 
such known bumps around the 6 + mass range; however, it is difficult to 
describe the background under a 9 signal due to these £ bumps. 

K s particles were reconstructed from secondary- vertex CTD tracks with 
transverse momenta pt > 0.15 GeV and pseudo-rapidities \r)\ < 1.75. 
The K° transverse momenta and pseudorapidities were required to have 
p T (K°) > 0.3 GeV and \r](K°)\ < 1.5. A very clean K° -> tt+tt" signal 
was obtained 2 . After requiring 0.483 < M(7r+7r~) < 0.513 GcV, the 
number of Kg candidates was ~ 867, 000 with only ss 6% background. 

Protons and antiprotons were selected from a wide dE/dx proton band, 
motivated by the Bethe-Bloch equation, defined for primary- vertex tracks 2 . 
Pion and kaon contamination was reduced by requiring the proton momen- 
tum to be less than 1.5 GeV and dE/dx to be above 1.15 minimum ionising 
particles (mips). The purity of the proton sample obtained was « 60%. 

The K s p(p) mass spectrum is shown in Fig. l(a-f) with a minimum Q 2 
ranging from 1 to 50 GeV 2 , as well as with Q 2 > 1 GeV 2 for two separate 
bins of the photon-proton centre-of-mass energy, W. For Q 2 > 10 GeV 2 
or Q 2 > 1 GcV 2 and W < 125 GeV, a peak is seen near 1.52 GcV. The 
histograms are the ARIADNE Monte Carlo (MC) simulation, normalised 
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Figure 1. M(K° s p(p)) for (a-d) Q 2 > 1, 10, 30, 50 GeV 2 ; (e-f) W < 125 and > 125 GeV; 
(g) Q 2 > 20 GeV 2 . The histograms are MC predictions normalised to the data above 
1.65 GeV. The solid line in (g) is a fit to the data using a background function (dotted 
line) plus two Gaussians (dashed lines). The inset shows the Kgp (open circles) and 
Kgp (black dots) combinations, compared to the combined sample fit scaled by 0.5. 

In Fig. 1(g) the Kg p (p) mass spectrum is shown for Q 2 > 20 GeV 2 
together with a fit to two Gaussians and a background of the form 
Pi(M - m p - m K a) P2 ■ (1 + P 3 (M - m p - m K o)), where M is the K%p 
mass, m p (m^o) is the proton (K°) mass and Pi, 2, 3 are free parameters. 
The fit x 2 /ndf (35/44) is significantly better than a one-Gaussian fit for 
the O only. The improvement is mainly in the low mass region, where the 
second resonance may correspond to the £(1480). The O peak position is 
M = 1521.5 ± 1.5(stat.)tl®(syst.) MeV, with a Gaussian width a = 6.1 ± 
1.6(stat.)^l'"(syst.) MeV, which is above, but consistent with the resolution 
(w 2 MeV). The fit gives 221 ±48 events above background, corresponding 
to 4.6<t. The probability of a fluctuation leading to the observed signal in 
the mass range 1.5 — 1.56 GeV is below 6 • 10~ 5 . Fitting the 9 with a Breit- 
Wigner convoluted with a Gaussian fixed to the experimental resolution, the 
intrinsic full width of the signal is estimated to be T = 8 ± A(stat.) MeV. 
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The signal is seen for both proton charges (inset in Fig. lg). The fitted 
number of events in the K^p channel is 96 ± 34. If the signal originates 
from the 0, this is a first evidence for the anti-pentaquark O - . 

The 6 production cross section was measured 4 in the kinematic region 
Q 2 > 20 GcV 2 , 0.04 < y < 0.95, p T (&) > 0.5 GeV and |ry(6)| < 1.5 to 
be a{ep -> eQ ± X -> eK°p ± X) = 125 ± 27 (stat.)±H(syst.)pb, where y is 
the lepton inelasticity. The acceptance was calculated using the RAPGAP 
MC, where S + baryons were treated as 6+ with M = 1.522 GeV, forced 
to decay 100% to Kgp(p). The G visible acceptance was « 4%. Fig. 2(a) 
shows the cross section integrated above several Q 2 min values. The cross- 
section ratio to that of A(1116), R = <r(9 + — > K p)/cr(X) (antiparticles 
are included), was measured in the same kinematic region. A baryons 
were measured in the decay mode A — > pn~ and protons were selected 
by dE/dx with identical cuts as for the 9 + . The A acceptance (« 10%) 
was calculated using the ARIADNE MC. The result for Q 2 > 20 GeV 2 is 
R = (4.2 ± 0.9±o g)%. Fig. 2(b) shows R for these Q 2 min values. It is not 
compatible with upper limits from HERA-B and ALEPH, where R < 0.5%. 
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Figure 2. (a) Visible cross sections for the Q ± baryon decaying to K°p(p) as a function 
of (b) Cross section ratio _R = ct(0+ — * K°p)/cr(A) as a function of Q^ in - 



If the observed 1.52 GeV peak is due to a / = 1 state, a ++ signal is 
expected in the K + p and K~p spectrum. Selecting protons and charged 
kaons using dE/dx, no peak was seen in the above distribution. A clean 10cr 
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A(1520) — > K p or K + p signal was seen with mass and width consistent 
with the PDG values 3 . The number of A(1520) and A(1520) are similar. 



3. Search for pentaquarks in the Sir channels 

The pp fixed-target NA49 Collaboration (y/s — 17.2 GeV) reported 5 obser- 
vation of the 2 multiplet pentaquark candidates and S°y 2 predicted 1 
at the bottom of the anti-dccuplet of baryons. They found narrow peaks 
in these Stt combinations at M rs 1862 MeV with a width < 18 McV. The 
significance of the signal for the sum of all 4 Bit channels is 5.8c. 

ZEUS searched for such states in its DIS HERA-I data 6 . 3 - (S+) 
states were reconstructed via the A7r~(A7r + ) decay channel, with 
A — > pir~ (A — * pir + ). Very clean A and S signals were obtained with 
130000 A + A and w 2600 5 + 3 candidates. In Fig. 3 the Stt invari- 
ant mass spectrum for Q 2 > 1 GeV 2 is shown. The left histograms show 
each charge combination separately. The right histogram is the sum of all 
Ett combinations. A clean 3° (1530) — * Ett signal of « 4.8cr is seen in the 
combined plot. No evidence is seen for the NA49 signal around 1862 MeV 
in any of the Ett mass plots. No signal is visible also with Q 2 > 20 GeV 2 . 
The discrepancy may be due to the fact that the ZEUS results come from 
the central rapidity region, while NA49 also covers the forward region. 
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Figure 3. M(B.tt) for Q 2 > 1 GeV 2 for each charge combination (left) and for all charge 
combinations combined (right). 



4. Search for a charmed pentaquark decaying to D*^p^ 

The existence of the strange pentaquark + implies that charmed pen- 
taquarks, — uuddc, should also exist. One type of model 7 predicts 
M(8°) 2710 MeV, which is below the threshold to decay strongly 
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to D mesons. Another model 8 predicts a which decays mainly to 
D~p or D°n (charge conjugate included) with M(0 9 ) = 2985 MeV and 
r(9°) w 21 MeV. If M(9°) is above the sum of the D* and p masses 
(2948 MeV), it can decay also to D^p*. 

The HI Collaboration found 9 a narrow signal in the D* ± p =F invariant 
mass at 3.1 GeV with a width consistent with the detector resolution. The 
signal was seen in a DIS sample of w 3400 D*± -> D°tt ± -> (K^tt ± )'k ± 
with a rate of w 1% of the visible D* production. A less clean signal of a 
comparable rate was seen also in the HI PHP sample. 

The 9 search of ZEUS in the D* ± p T mode was performed with 
the full HERAT data 10 . Clean D* ± signals were seen in the 
AM = M(D* ± ) - M(D°) plots (Fig. 4 left). Two D*± -» D ^ 
decay channels were used with D° — » X =F 7r ± and I? — > if T 7r ± 7r + 7r _ . 
The 0° search was performed in the kinematic range \ri(D*)\ < 1.6 
and p T (D*) > 1.35 (2.8) GeV and with AM values between 
0.144 - 0.147(0.1445 - 0.1465) GeV for the Ktttt (Ktttttttt) channel. In 
these shaded bands a total of w 62000 D* } s was obtained (Fig. 4a-b left) 
after subtracting wrong-charge combinations with charge ±2 for the D° 
candidate. Selecting DIS events with Q 2 > 1 GeV 2 yielded smaller, but 
cleaner D* signals with a total of w 13500 D*'s (Fig. 4c-d left). 

Protons were selected with pr(p) > 0.15 GeV. To reduce the pion and 
kaon background, a parameterisation of the expected dE/dx as a function 
of P I m was obtained using tagged protons from A decays and tagged pions 
from Kg decays. The \ 2 probability of the proton hypothesis was required 
to be above 0.15. Fig. 5 shows the M(D*p) = M(Kmrp) ~ M(Ktttt) + 
M(D*)pdg distributions for the Kirn channel for the full (left) and the 
DIS (right) samples, where M(D*)pdg is the D* ± mass 3 . In the low- 
P selection (Fig. 5b), a clean proton sample separated from the 7r and 
K dE/dx bands was obtained by taking only tracks with P < 1.35 GeV 
and dE/dx > 1.3 mips. In the high-P selection (Fig. 5c) only tracks with 
P{p) > 2 GeV were used. The latter selection was prompted by the HI 
observation 9 of a better 6° signal-to-background ratio for high proton 
momenta. No narrow signal is seen in the Ktttt (Fig. 5) as well as in the 
Kmnnr (Fig. 4b, d right) channel. The Ktttt analysis was repeated using 
very similar selection criteria as in the HI analysis 9 . No indication of a 
narrow resonance was found in either the DIS or the PHP event sample 10 . 

95% C.L. upper limits on the fraction of D* mesons originating from 
6° decays, i?(9 9 — » D*p/D*), were calculated in a signal window 
3.07 < M(D*p) < 3.13 GeV for the Ktttt and K tttt-ktt channels. A visible 
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Figure 4. Left: AM distributions (dots) for (a) D* — » K-ktt and (b) D* — » Kttttitit 
candidates. Events with Q 2 > 1 GeV 2 for the two channels, respectively, are shown 
in (c) and (d). The histograms are for wrong charge combinations. Right: M(D* ± p = F) 
distributions (dots) for the same samples. Solid curves are fits to a background function 
(see text). Shaded historgams are MC 0^ signals, normalised to 0j?/D* = 1%, on top 
of the background fit. 



rate of 1% for this fraction (Fig. 4 right), as claimed by HI 9 , is excluded 
by 9cr (5a) for the full (DIS) combined sample. The M(D*p) distributions 
were fitted to the form x a e~ bx+cx , where x = M(D*p) — M(D*) - m p 
(Fig. 4 right). The number of reconstructed 0° baryons was estimated by 
subtracting in the signal window the background function from the observed 
number of events, yielding R(Q° C -> D*p/D*) < 0.23% and < 0.35% for the 
full and DIS combined two channels. The acceptance-corrected rates are, 
respectively, 0.37% and 0.51%. The 95% C.L. upper limit on the fraction 
of charm quarks fragmenting to 0° times the branching ratio 0° — ■> D*p 
for the combined two channels is /(c — * 0^) -B@o_, D , p < 0.16% (< 0.19%) 
for the full (DIS) sample. 

5. Summary 

The ZEUS HERA-I data sample was used to search for narrow baryonic 
pentaquark candidates. For the inclusive DIS sample a 4.6cr narrow signal 
was seen in the fragmentation region in the combined M(Kgp) and M(Kgp) 
plot at the mass range. If due to the baryon, this is the first evidence 
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Figure 5. Left: M(D* ± p = F) distributions for the i^7T7r channel (dots) with (a) all proton 
candidates, (b) candidates with P(p) < 1.35 GeV and dE/dx > 1.3, and (c) candidates 
with P(p) > 2 GeV. Histograms show the M(i?* ± p ± ) like-sign combinations. Right: 
Same for DIS events with Q 2 > 1 GeV 2 . 



for the anti-pentaquark 6 . The cross-section ratio a(Q + — > K°p)/a(A) 
for Q 2 > 20 GeV 2 is (4.2 ± 0.9±H)%. No evidence is found for the NA49 
S-7T signal at 1862 MeV in the inclusive DIS sample. No resonance structure 
is seen in MiD^pT) around 3.1 GeV. The 95% C.L. upper limit on the 
visible rate R{Q° C D*p/D*) is 0.23% (0.35% for DIS). The ZEUS data 
are not compatible with the HI result of « 1% of the above rate. Such a 
rate is excluded by 9a for the full data and by 5cr for the ZEUS DIS data. 
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